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A simple approach is presented for the modeling of complex oscillatory thermal-fluid systems capable of 
converting low grade heat into useful work. This approach is applied to the NIFTE, a novel low temperature 
difference heat utilization technology currently under development. Starting from a first-order linear 
dynamic model of the NIFTE that consists of a network of interconnected spatially lumped components, the 
effects of various device parameters (geometric and other) on the thermodynamic efficiencies of the device 
are investigated parametrically. Critical components are highlighted that require careful design for the 
optimization of the device, namely the feedback valve, the power cylinder, the adiabatic volume and the 
thermal resistance in the heat exchangers. An efficient NIFTE design would feature a lower feedback valve 
resistance, with a shorter connection length and larger connection diameter; a smaller diameter but taller 
power cylinder; a larger (time-mean) combined vapor volume at the top part of the device; as well as 
improved heat transfer behavior (i.e. reduced thermal resistance) in the hot and cold heat exchanger 
blocks. These modifications have the potential of increasing the relevant form of the second law efficiency 
of the device by 50% points, corresponding to a 3.8% point increase in thermal efficiency. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

With ever-increasing environmental concerns, including that of 
climate change, but also of energy security in the light of finite 
resources of common fossil fuels, it is becoming increasingly impor¬ 
tant to consider alternative clean and sustainable energy solutions. 
The efficiency and power density of appropriate thermodynamic 
systems are key quantities on which their usefulness and wider 
application rely. This is especially true for devices that are capable of 
operating with marginal energy sources and close to ambient 
temperatures. In this paper we present a modeling framework based 
on thermodynamic, heat transfer and fluid-mechanical principles for 
the early stage development of oscillatory low grade heat (e.g. solar 
energy or waste heat) utilization systems. We demonstrate the 
approach by applying it to a promising new low grade heat utilization 
technology that has attracted attention recently as a result of 
potential reliability advantages as well as reduced capital and oper¬ 
ating costs, on account of its few moving parts relative to mechanical 
heat engines. Specifically, we consider the effects of various system 
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parameters on a number of key efficiency definitions that can be 
formulated in order to describe the performance of the system. 

The “Non-Inertive-Feedback Thermofluidic Engine” (NIFTE), as 
proposed in Refs. [1—3], can be described as a “two-phase 
unsteady heat engine” in which persistent and reliable thermo¬ 
dynamic (pressure, temperature, etc.) oscillations are generated 
and sustained by steady external temperature differences. These 
desirable oscillations are driven by and give rise to heat and fluid 
flows, which involve the evaporation (boiling) and condensation 
of the working fluid. Hence, the NIFTE can also be considered 
a two-phase realization of a class of devices known as “thermo¬ 
fluidic oscillators”, which includes thermoacoustic engines [4—7], 
dry free-liquid-piston (Fluidyne) engines [8—11], free-piston Stir¬ 
ling engines [12—14], pulsejets and pulse-tubes [15—19]. In 
common with other two-phase thermofluidic oscillators, the 
NIFTE is particularly well suited to the conversion of low grade 
heat. It can be used to produce useful (hydraulic) work for fluid 
pumping, heating and/or cooling and niche power generation 
applications; in a variety of solar thermal and energy recovery 
settings, and in combined heat and power-pumping (CHP) 
schemes. The NIFTE has been demonstrated as being capable of 
operating across temperature differences down to 30 1< between 
the heat source and sink. 
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Nomenclature 

C Electrical capacitance equivalent (m 4 s 2 kg -1 ) 

E (—* P) Electrical potential/voltage equivalent 
(Pa) 

I (—>U) Electrical current equivalent (m 3 s _1 ) 

L Electrical inductance equivalent (kg itT 4 ) 

P Pressure (Pa) 

Q Heat flow (J) 

R Electrical resistance (kg m 4 s 1 ) 

S Entropy (J K 1 ) 

T Temperature (°C or I<) 

U = V Volumetric flow rate (m 3 s -1 ) 

V Volume (m 3 ) 

W Work done by a thermodynamic system 

(J) 

X Exergy 0) 

Z Electrical impedance equivalent (kg ur 4 s _1 ) 

c p Specific isobaric heat capacity at constant pressure 

J kg ' K ') 

c v Specific isochoric heat capacity at constant volume 

(J kg- 1 K- 1 ) 
d Diameter (m) 

/= 1/t Frequency of a cycle or device (Hz) 
g = 9.81 Gravitational acceleration (m s~ 2 ) 
h Enthalpy per unit mass of working fluid, or specific 

enthalpy (J kg- 1 ) 

k Thermal conductance or inverse of the thermal 

resistance (W I< “') 

1 Length (m) 

m Mass (kg) 

s Entropy per unit mass of working fluid, or specific 

enthalpy (J kg -1 I< -1 ) 

u Bulk flow velocity averaged over flow cross-section 

(m s- 1 ) 

v = \/p Volume per unit mass of working fluid, or specific 
volume (m 3 kg- 1 ) 

Creek 

g Dynamic or absolute viscosity (kg nrr 1 s _1 ) 

p Density (kg m -3 ) 

i=l// Cycle period (s) 

to = 2t if Angular frequency (rad s _1 ) 


Dimensionless 

G Linear system gain (-) 

Re Reynolds number (-) 

Wo Womersley parameter (-) 

c d Dimensionless friction factor/drag coefficient ( -) 

y = c p /c v Ratio of specific heat capacities (-) 

t] Thermodynamic efficiency (-) 

Operations 

<(•)> 

(0 

(•) = 

(•)' = (•)-(■) 

( ; )= d(-)/dt 
(•)* 

*{-.•} 

A(-) 

Other 

ad Adiabatic process or adiabatic vapor chamber 

ao Mean value during an adiabatic process in the 

adiabatic vapor chamber 
c Cold, or heat sink 

d Displacer cylinder 

dev Device, or overall/global to the complete fluid 

pumping device 
ex Exergetic 

f Feedback, or feedback connection and valve 

fg Change over an evaporative phase change 

process 

h Hot, or heat source 

in Into the thermodynamic system over a cycle 

1 Load 

o Mean/time-averaged value 

out Out of the thermodynamic system over a cycle 

p Power cylinder 

s Common designation for heat source and heat 

sink 

sat For a fluid in the saturation region 

th Thermal, or thermal domain involving the 

heat exchangers 

wfl Working fluid, or specific to the working fluid 

thermodynamic cycle 


Averaged, as stated 
Time-mean average 
Laplace transform of 

Temporal fluctuation of a quantity around its 
mean 

Rate of change, or flow rate of 

Value of a parameter relative to its nominal or 

maximum investigated value 

Phase angle between the stated variables 

Change/difference in 


Thermofluidic oscillators have many dynamic similarities with 
analog electronic oscillator circuits, and thus, electrical analogies 
have been used to predict approximate stability/instability criteria 
and to estimate first-order heat and work flows, and efficiencies. 
The analogies, used extensively by Backhaus and Swift [4] in the 
analysis and development of thermoacoustic engines, were 
extended to include a description of exergy (or availability) flows 
and to allow for exergy losses due to irreversible heat transfer in 
order to model the N1FTE by Smith [lj. In this work a “conven¬ 
tional” exergetic efficiency [20] was introduced as the second law 
efficiency most appropriate to thermodynamic models described 
by electrical analogies in which the heat source and heat sink 
temperatures are represented as an externally imposed linear 
temperature gradient rather than separate constants. This exergetic 
efficiency, defined below in Section 2.4, is independent of absolute 
temperatures as well as temperature differences (and amplitudes). 
It is therefore conveniently invariant to changes in operating 


conditions and was proposed as a general figure-of-merit for the 
characterization of the thermodynamic reversibility of the engine, 
and by extension, its performance. 

The electrical analogy model of Smith [1—3] showed good 
agreement with experimental trends obtained from measurements 
on an actual NIFTE prototype over a certain range of conditions and 
prototype design parameters, though it also showed non-negligible 
deviations in some cases [2,3]. Considerable effort is placed in these 
references to explain the underlying reasons for these deviations, 
including for example the existence of shuttle (irreversible heat 
transfer) loss mechanisms in non-heat exchanging (or, nominally 
adiabatic) components, and the presence of fluid (liquid) inertia in 
the real NIFTE. Unfortunately, the original model was not extended 
to include such effects, and in particular, although mentioned in 
this early work the effect of inertia was not studied. 

Yet, any real embodiment of such a system will encompass 
inertial effects to some extent due to the finite mass of any 
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contained liquids. A dependence on inertia in order to allow the 
system sustained pressure oscillations with large amplitudes, 
which is vital for the achievement of high power densities, implies 
a requirement for bulky tuning lines and resonators that the 
designer would like to avoid. An understanding of the need for and 
of the effects of finite inertia is crucial for the design of the next 
generation of these novel and complex systems, both in terms of 
enhanced efficiency and performance. The effects of inertia were 
incorporated into the original “non-inertive” electrical analogy 
model of Smith [1], by Solanki et al. [21]. The current work uses this 
latter “inertive” model as a starting point. 

In Solanki et al. [21 ] the authors studied the effect of introducing 
(and parametrically varying the degree of) inertia on key NIFTE 
performance indicators, such as the minimum temperature 
gradient between the heat source and sink necessary for operation 
(continuous oscillations) and the resulting oscillation frequency. 
The revised inertive model resulted in a much improved prediction 
of these indicators relative to the original non-inertive equivalent, 
over a wider range of values. This suggests that the inclusion of 
inertia is indeed essential for the correct capture of the observed 
behavior of the NIFTE. The inertive model was then used to 
investigate the expected efficiency of the NIFTE for a typical and 
realistic set of system properties and design parameters. 

However, this previous study focused only on the exergetic 
(second law) efficiency of the engine and did not present results for 
its thermal (first law) efficiency. In addition we note that, in general, 
the real NIFTE is an energy producing device that can have a different 
(indeed, reduced) global efficiency from the efficiency of its enclosed 
thermodynamic cycle, owing to processes of parasitic dissipation 
that are external to the working fluid cycle, but still internal to the 
overall device. In fact, the NIFTE depends on a throttle valve (referred 
to as the “feedback valve”) to create a phase shift between the heat 
flow and the power stroke. The existence of this phase shift results in 
a parasitic fluid-mechanical viscous loss that leads to a difference 
between the efficiency of the thermodynamic cycle undergone by 
the working fluid, and the overall or global efficiency of the complete 
device. In Smith [1—3] the effect of this component on the exergetic 
efficiency was assumed negligible. In Solanki et al. [20], it was 
included throughout. However, only its effect on the overall exer¬ 
getic efficiency of the whole engine was considered. Whilst this 
provides some useful insights, increases in exergetic efficiency may 
lead to decreases in overall thermal efficiency, since decreases in 
thermal dissipation and parasitic losses are often accompanied by 
a reduction in the temperature difference between heat addition 
and heat rejection. Hence, a more comprehensive approach is 
required for a complete understanding of the NIFTE. 

In this paper we use the extended inertive model of the NIFTE 
[21 ] to investigate the effect of the fluid-mechanical parasitic losses 
both on the exergetic efficiencies and on the thermal efficiencies of 
the NIFTE pump. We refer to efficiencies in which the parasitic 
fluid-mechanical losses are ignored as “working fluid” efficiencies, 
and efficiencies in which they are included as “device” efficiencies. 
Therefore, with regards to the NIFTE, one needs to distinguish 
between four measures of efficiency: thermal and exergetic overall/ 
global device efficiencies; and thermal and exergetic working fluid 
efficiencies. Solanki et al. [21] restricted their attention to the 
exergetic overall device efficiency, while in Refs. [1—3] the parasitic 
losses were assumed to be negligible altogether (i.e. no distinction 
was made between device and system efficiencies). Here, we 
investigate the consequences of this assumption and show that it is 
not valid in certain regions of the parameter space. We use the 
results of this process to construct a framework for including 
parasitic losses in future thermofluidic oscillator models, in order to 
gain a more complete understanding of the thermodynamic 
performance potential of the NIFTE. 


2. Methods 

2.3. Problem formulation: heat, power, exergy and efficiency definitions 

Consider an arbitrary, externally heated work-producing device, 
within which power is generated by a working fluid that undergoes 
a thermodynamic cycle with frequency/. Rather than dealing with 
thermodynamic states we will deal with flow rates of thermody¬ 
namic properties (e.g. “entropy flow rate” S), so that irreversible 
processes may be quantified legitimately as increases in entropy 
flow rate. Lower case letters are used to indicate specific quantities 
(per unit mass of working fluid), such that for instance s is the 
specific entropy flow rate. 

The temperature of the working fluid of the device that 
undergoes this “main cycle” can be written as a decomposition 
T=T 0 + T, where T are fluctuations about the time-averaged value 
T 0 . By definition the time-average of the fluctuations is zero T = 0. 
The working fluid is heated by an external source and cooled by an 
external sink, through contact with appropriate heat exchangers. 
The source and sink on the other side of the heat exchangers 
undergo separate processes. Let us denote by Th the temperature of 
the source from which the main cycle obtains heat, and by T c the 
temperature of the sink to which it rejects heat. Assuming that the 
main cycle is heated and cooled by exclusively irreversible heat 
transfer, i.e. pure Fourier conduction (diffusion) and Newtonian 
convection (surface) phenomena, Th and T c are bounded by lines of 
constant heat flow rate Q_ that set a limit to the extent of the main 
cycle, wherein T = Q/S. In addition, let AS = f Sdt > 0 be the 

s>o 

entropy change during the process of heat addition (S > 0), such 
that the entropy change during heat rejection (S < 0) is equal to 
this in magnitude but opposite in sign (AS < 0); and let V be the 
volume of the working fluid, such that LI = V is the volumetric 
displacement flow rate responsible for power generation. 

Then, based on an arbitrary datum state (e.g. chose the time- 
averaged T 0 , P 0 ). the rate of heat input to the main cycle Q in . the 
net power produced W, and the “net exergy (or equivalently 
availability) flow rate” X are: 

Q ln = J T(S)dS = fT 0 \S+ J T'(S)dS (1) 

S,AS>0 S,AS>0 

W = <f> P(V)dV = j P(U)dU = Q = Qin - Q out = j T(S)dS 

( 2 ) 

X ln = J [T h (S) -T 0 ]dS; X ou t = J [T C (S) -T 0 dS]; 

S,AS>0 S,AS<0 

x = x in -Xout = ^r s (s)ds (3) 

where T s is used to denote either Th or T c during the sequential 
heating (S >0; AS > 0) and cooling (S <0; AS < 0) parts of the 
cycle. In addition, we define the equivalent rate of heat gained from 
the heat source Q h external to the main cycle as, 

Qh = J T h (s)ds = fT 0 \s+ J r h (s)ds. (4) 

S,AS>0 S,AS>0 

which from Eq. (3) is actually the thermal exergy input X in plus the 
constant (time-mean component) /T 0 AS supplied by the heat 
source during the heat addition stage. 

From Eq. (2) we can see that the power produced by the device 
VV = Q may be described by the area inside a T-S diagram that 
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represents the internally reversible part of the main thermodynamic 
cycle, multiplied by the frequency of the cycle/. This area is enclosed 
by a larger area that, when multiplied by frequency/ represents the 
net exergy flow rate X of the full irreversible cycle, with the differ¬ 
ence between these areas representing power loss through irre¬ 
versible processes. Hence, X is the maximum power available to the 
cycle, given the external boundary conditions. It is clear from Eqs. (2) 
and (3) that changes in T 0 will have no effect on W orX. Conversely, 
from Eqs. (1) and (4) we observe that both Q in and Q h do indeed 
depend on T 0 . and consequently on the magnitude of relative to T 0 . 

Next we consider two classes of efficiency: (1) the thermal 
efficiency ?)th, defined as r] th = VV/Q h ; and (2) the exergetic effi¬ 
ciency rj ex , defined as y ex = W/X. In earlier models of the NIFTE the 
parasitic fluid-mechanical viscous dissipation in the feedback valve 
was assumed negligible [1—3], We refer to efficiencies based on this 
assumption as (thermal or exergetic) “working fluid” efficiencies 
?j w fi and efficiencies for which this fluid-mechanical dissipation is 
accounted for as (thermal or exergetic) “device” efficiencies Tjdev- In 
this paper we consider the effect of various device parameters on 
these efficiencies, with reference to the NIFTE which is introduced 
in more detail below. 

2.2. Non-Inertive-Feedback Thermofluidic Engine (NIFTE) 

A schematic of the NIFTE can be seen in Fig. 1. Details of the 
construction and operation of the NIFTE can be found elsewhere 
[3,21], but briefly the NIFTE comprises two vertical cylinders (1,2), 
connected at the top (3) and bottom (4) with horizontal pipes. The 
whole interconnected volume contained within these main 
components (chambers and tubes) is filled with a single working 
fluid that exists in both the vapor (white) and the liquid (gray) 
phase at the same time. The liquid and vapor volumes are separated 
by two vapor—liquid interfaces denoted in Fig. 1 by levels (6) and 
(7), the first of which is inside the power cylinder on the left (1) and 



Fig. 1. Schematic of the NIFTE. The area in white corresponds to the vapor volume in 
which the working fluid is in the vapor phase, while the area in gray corresponds to the 
liquid volume in which the working fluid is in the liquid phase. Component 1 is the 
power cylinder (denoted by subscript 'p' in the text) and 2 is the displacer cylinder 
('d'), 3 is the vapor connection, and, 4 and 5 are the liquid feedback connection and 
valve ('f). Levels 6 and 7 are the working fluid vapor-liquid interfaces in the power 
(left) and displacer (right) cylinders. The combined vapor volume 8 (above levels 6 and 
7) is assumed to be an adiabatic vapor chamber {‘ad’). Interface 9 is connected to the 
load line (T). Components 10 (hot) and 11 (cold) are heat exchanger blocks (‘th’). 


the second is inside the displacer cylinder on the right (2). The total 
volume of vapor at the top of the NIFTE (8) includes the vapor 
volumes in the power and displacer cylinders, together with the 
volume of the horizontal “vapor connection” pipe (3). This 
combined space is assumed to be an “adiabatic vapor chamber” (8). 
The second horizontal pipe that connects the power and displacer 
cylinders further down is referred to as the “liquid feedback 
connection” (4) and contains a valve termed the “feedback valve” 
(5). The power cylinder (1) transmits the fluid displacement to the 
load (9) below it. We treat the case where the working fluid and 
power transmitting fluid are the same. 

The displacer cylinder (2) contains two heat exchanger blocks 
arranged vertically, namely the hot (10) and cold (11) heat 
exchangers. The heat transfer coefficient associated with single phase 
convection is much lower (typically by orders of magnitude) than that 
associated with phase change, and hence we neglect all modes of heat 
transfer other than evaporation and condensation. Evaporation 
occurs when working fluid in the liquid phase comes into contact with 
the surface of the hot heat exchanger (10), whereas condensation 
occurs when working fluid in the vapor phase comes into contact with 
the cold heat exchanger (11). The latter is the dominant heat transfer 
process at the instant portrayed in Fig. 1, when the vapor—liquid 
interface (7) is at the height of the cold heat exchanger (11). 

Starting from the instant of condensation depicted in the figure, 
a reduction of mass of working fluid in the vapor phase (white) leads 
to a decrease in pressure in the vapor volume of the NIFTE (i.e. the 
combined space at the top of the power and displacer cylinders (8), 
and in the vapor connection). This reduction in pressure causes level 6 
to be forced upwards in the power cylinder (1), and fluid to be sucked 
upwards into the NIFTE from the load line (9). As level 6 rises in the 
power cylinder (1), the hydrostatic pressure difference between the 
power (1) and displacer (2) cylinders first decelerates and then 
reverses the upwards motion of level 6. At this crucial point in the 
operation of the NIFTE the feedback connection (4) allows liquid to 
flow from the power (1) to the displacer cylinder (2) that contains the 
heat exchanger blocks. This, in conjunction with the reduced pressure 
in the vapor phase (8), causes level 7 to rise, bringing liquid into 
contact with the hot heat exchanger (10). The evaporation of this 
liquid increases of volume of vapor (white) and leads to pressuriza¬ 
tion of the vapor volume that forces level 6 downwards, while fluid is 
displaced downwards and out of the NIFTE into the load line (9). As 
level 6 falls, the hydrostatic pressure difference between the power 
(1) and displacer (2) cylinders first decelerates and then reverses the 
downwards motion of level 6. The feedback connection (4) now 
allows liquid to flow from the displacer (2) to the power (1), while the 
pressurization causes level 7 to descend onto the cold heat exchanger 
whence condensation can re-commence, completing one cycle of the 
NIFTE. The process then repeats itself. The role of the feedback valve 
(5) is, importantly, to regulate the flow rate through the feedback 
connection (4), which determines the rate of liquid exchange 
between the power (1) and displacer (2) cylinders, and consequently 
the oscillation frequency of the entire device. The oscillating fluid 
displacement into and out of the load line (9) can be harnessed in 
various ways, one of which involves the use of a pair of check valves in 
opposite directions and in parallel to each other in order to achieve an 
uni-directional flow for pumping [3]. A linear approximation of the 
operation of the NIFTE can be described by an LRC circuit, that is, an 
electrical circuit with inductors (I), resistors ( R ) and capacitors (C), as 
explained in the following section. 

2.3. Spatially lumped and one-dimensional linear dynamic model 

In Smith [1—3] a simple, but powerful model for the dynamic 
behavior of the NIFTE was proposed, following Backhaus and Swift 
[4], Ceperley [7] and Huang and Chuang [15]. The model of the 
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Pi 


Fig. 2. N1FTE model circuit diagram, where C is the feedback gain, R, is a resistance, C, is a capacitance, I, is an inductance, P, is a pressure and (J, is a volumetric flow rate. Subscript 
‘th’ denotes the thermal domain (heat exchanger blocks 10 and 11 in the repeat of Fig. 1 on the left), ‘ad’ the adiabatic domain (vapor volume 8, above liquid levels 6 and 7), T the 
load (below 9), 'p' and 'd' the power (1) and displacer (2) cylinders, and 'f the feedback connection: pipe and valve (4,5). 


NIFTE involved a modification to the earlier models by these 
authors to incorporate heat flows with irreversible heat transfer, 
and an extension that allowed a description of phase change 
processes. A description of this model is presented in Ref. [3], 
however, the essential elements are reproduced in a more 
condensed form below, with some additional insights. 

The model involved lumped (spatially averaged and thus inde¬ 
pendent) and first-order linearized sub-models for each compo¬ 
nent of the NIFTE that were interconnected to derive a complete 
model for the whole device. The linearization allows analogies to be 
drawn with analog electrical components, thus enabling an elec¬ 
trical network to be constructed for the NIFTE. The network 
contains resistors (accounting for viscosity, fluid drag and thermal 
resistance) and capacitors (for gravity and compressibility). Here 
we extend this approach, with the added inclusion of inductors to 
capture finite inertial effects in line with Solanki et al. [21], 

In more detail, the dominant physical process undergone in 
each sub-component is firstly identified and modeled linearly. By 
drawing analogies between the physical variables of pressure (P) 
and voltage (£), volumetric flow rate (Lf) and current (/), tempera¬ 
ture (T) and voltage (£), and entropy flow rate (S) and current (/), the 
linear, lumped equations involving the raw thermal-fluid variables 
are transformed into a suitable network of electrical elements. The 
use of voltage to represent both pressure and temperature, and the 
use of current to represent both volumetric flow rate and entropy 
flow rate necessitates the definition of coupling equations, if both 
are to be used within a single model. 

Referring to the electrical network in Fig. 2, resistances to heat 
or fluid flow are represented by a resistors (R); hydrostatic pressure 
(i.e. gravitational potential energy) of the liquid and adiabatic 
compressibility of the vapor are both represented by capacitors (C); 
and fluid inertia is represented by inductors (I). The dynamic 
equation of each element is / = £/R, I = C.E and / = E/I, respec¬ 
tively, where £ is the potential difference across the component and 
/ is the current through it. The physical quantities used to evaluate 
the value of each electrical component in Fig. 2 are shown in 
Table 1, with an explanation of all variables in the caption. 


Note that the network in Fig. 2 represents a model of connected 
interacting processes, and that some of these are thermal (heat) 
flow processes while others are fluid flow processes. Now, the rate 
of change of pressure can be written as, 


dPad 
d t 


■Sfh 



-U ad 

Vad 



•Sth 


(5) 


The first term is determined purely by heat transfer processes and 
the second by adiabatic ones. For a two-phase process, it can be 
shown without approximation from the Clausius—Clapeyron 
equation, and by noting that PlA is constant for a reversible adia¬ 
batic (isentropic) process, that this is equivalent to, 


dPad 

df 


‘fg 


(Cv)v fg 


Sth + 


yPad 

Vad 


U ad 


( 6 ) 


where (c„) is the mass-averaged specific heat capacity at constant 
volume of a volume of liquid—vapor mixture in equilibrium. With 
reference to the second term on the right-hand side (RHS) of Eq. (6), 
consider the adiabatic compression/expansion process undergone 
by the working fluid vapor in the vapor chamber (8 in Fig. 1) due to 
the time-varying changes to the liquid levels in the power (6) and 
displacer (7) cylinders. This is a ‘fluid’ process and as such can be 
purely described by changes in a pressure (P) and a volumetric flow 
rate (LI). In this case, pressure is the ‘potential’ (or ‘voltage’) and the 
volumetric flow rate is the ‘current’. Let P a d, V a d be the time-varying 
pressure and volume of the vapor, such that the process can be 
described by P ad Vj d = cons,: - F° r sma H perturbations of pressure 
and volume around an equilibrium P ao = P 0 . Vao we have, 


1 dP ad y dU ad 
Pao df k/o dt 


=>u ad (t) 


^ao dP ad (t) 
yPao dt 


kao 
TPai 


sPad(S), 


Vao 

yP ai 


Pad (f) °r U id {S) 


(7) 


Table 1 

Expressions for the components appearing in the NIFTE as modeled in Fig. 2. The subscript ‘fg’ refers to phase change and ‘o’ to a time-averaged value; ASf g and Avf g denote the 
specific entropy and specific volume changes associated with phase change; k is a relevant thermal conductance; T 0 and P G = P ao are the time-averaged temperature and 
pressure; V ao is the mean total volume in the vapor phase; p, and p are the viscosity and density of the liquids; and /, and d, are pipe lengths and diameters. 


Thermal-Fluid Effect 

Electrical Component 

Analogy 

Expression 

Heat exchanger thermal resistance 

Resistor 

R th = ( s fg/ v fg) 2 Po/k 

Feedback valve & Load flow drag (viscous/pressure) 

Resistor 

R, = 128 /Ttdf &R f = 128/ilf/ndf 

Power & Displacer cylinder hydrostatic pressure 

Capacitor 

Cp = Ttdl/4pg&Ci = -Kd\j4pg 

Adiabatic vapor chamber compressibility 

Capacitor 

Qd = Vao/yPao 

Power & Displacer cylinder inertia (liquid mass) 

Inductor 

L p = 4plp/ndp & L d = 4pl d /Tidj 

Feedback connection & Load inertia (liquid mass) 

Inductor 

Lf = 4plf/Tcdf & Lf — 4pl\/ndf 
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that behaves equivalently to a capacitor / = C ad -£, with Qd = 
Vao/yPao- Note that here the tilde signifies the Laplace transform of 
the function, e.g. H ad = i?{H ad }, and s is the Laplace variable. 

In another example, not considered by Smith [1—3], we assume 
that the fluid flow in the load (9) is akin to viscous laminar flow in 
a smooth pipe (i.e. at low Reynolds numbers Re) and at low 
frequencies (i.e. at low Womersley parameters Wo) when the flow 
is quasi-steady. Let Ft be the pressure difference across the load and 
l/i the volumetric flow rate through it. In the absence of inertia 
(Lj = 0 in Fig. 2) Pi acts only across R h and the flow resistance (i.e. 
dimensionless friction factor/drag coefficient c d ) is modeled as, 


assume that T 0 s fg >> (c„)AT such that (c 1/ )^ + s fg =s fg (correct to 

* O 

about 97% for water/steam and 87% for pentane at typical N1FTE 
pressure amplitudes), we obtain, 

s th = ^f; th . (11) 

v fg 

The second coupling relation concerns the two potential variables. 
Continuing with the same assumption of thermal processes 
dominated by phase change used above for the flow variable, we 
have, 


r _ P i ; 1 _ 64 -.rj 1 p fs1 

d 1 2^1 Re d 1 128,u/i/7Td4 Pl ’ 

2 P] 

that resembles a resistor / = 1/Ri ■£, withRi = 128 (il^/nd^. Here, 
ill = Ui/(7td^/4) is the bulk velocity in the load. At higher 
frequencies one can no longer neglect inertia, which is then 
included in the form of an inductor Li in Fig. 2. This component is 
sufficient to account for changes in the amplitude and the phase of 
Pi relative to Hi, as demonstrated successfully in Ref. [22]. 

Now, consider instead the section indicated by the enclosed box 
in Fig. 2. This represents the heat exchangers (10 and 11), and 
specifically, the process of two-phase heat transfer by which the 
periodic lapping of the hot/cold heat exchanger blocks by the 
working fluid liquid level (7) in the displacer cylinder (2) drives an 
alternating process of generation/condensation of vapor directly 
into/out of the adiabatic vapor chamber (8). This is a ‘thermal’ 
process that can be described by changes in temperature (T) and 
entropy flow rate (S), whereby the ‘potential’ (or ‘voltage’) is now 
the temperature and the ‘current’ is the entropy flow rate. To first- 
order this process is described by, 

S|h ~T7/k lT ’~ Tl (9 > 

where k is the thermal conductance (inverse of the thermal resis¬ 
tance) between the solid block and the working fluid, and T 0 is the 
mean boiling point temperature over which the working fluid 
evaporates and condenses. 

It should be stated that for the purposes of our linearized model 
T 0 = T sat is the saturation temperature that corresponds to a satu¬ 
rated vapor pressure of the working fluid equal to the time- 
averaged pressure in the adiabatic vapor chamber 
Pad = Rao = Psat- It is also equal to the mean temperature between 
the hot source and cold sink, as well as the time-averaged working 
fluid temperature in the vapor phase within the adiabatic vapor 
chamber T 0 = T sat = T s = T. 

Returning to the section indicated by the enclosed box in Fig. 2, 
we note that the thermal process variables must be coupled to (or 
‘referred to') the fluid variables in the rest of the circuit. For this, we 
require a pair of coupling equations. The first one is for the flow 
variable, 



( 12 ) 


where (3T/3P) sat is the rate of change of temperature with pressure 
in the saturation region. We note further to Eq. (10) that after 
neglecting bulk kinetic and potential energy changes, energy 
conservation requires that, 


Q. — H=>TS t h — Pad^th + ^th (Ci/) AT 
=>S — Pad [/ (Cy) AT 

Sth - "T th + ^TT tlv 


(13) 


The first term on the RHS of Eq. (13) represents the displacement 
work advected away from the solid—fluid interface and the second 
term represents the internal energy advected from the interface. 
Eqs. (10) and (13) can be combined to give (3T/3P) sat = Vf g /sf g , which 
is correct to >99% for pentane at typical NIFTE pressure amplitudes. 

Having derived the coupling equations we return to Eq. (9), 
which describes the process of two-phase heat transfer in the heat 
exchangers. This process can now be re-written as, 


a r 

Nth = (Pth ~ Pad) 

^ To/k 

v fg 


. 2 (%-Pad), (14) 

^Sfg/tifgJ To/k 


and the behavior can be modeled by a resistor / = l/R th -£, with 
Pth = ( s fg/ ! 'fg) To/k. 

Results concerning all components in Fig. 2 are summarized in 
Table 1. Finally, it should be stated that the choice of all thermal- 
fluid-electrical analogies has been made to enforce the condition 
that electrical power, which can only be dissipated in resistors, is 
directly equivalent to either: (1) a rate of thermal exergy dissipated 
in the thermal domain due to irreversible heat transfer 
/•A£ = H AP = S AT = X (here, H and AP are the thermal quan¬ 
tities S and AT referred to the fluid domain); or (2) a rate of exergy 
dissipated in the fluid domain / A£ = H AP = W, either by 
conversion to useful work done on a load, or by conversion to lost 
work in a parasitic component. 


2.4. Solution procedure 


■Sth = AT + rh th s f g, (10) 

where m th the mass flow rate of fluid over the solid boundary that 
changes phase, (c„) is a suitably averaged specific heat capacity at 
constant volume of the liquid—vapor mixture, AT is the total 
temperature change of the liquid—vapor mixture including excur¬ 
sions outside the saturation region, and Sf g the specific entropy 
change during evaporation. Given the large density change asso¬ 
ciated with phase change, the net mass flow away from the solid 
surface can be re-expressed as m th = H th / v fg . For small oscillations 
about P ao , T 0 in the saturation region of the working fluid, and if we 


Similarly to Ref. [21], a “nominal” value of each electrical 
parameter in Fig. 2 was calculated in the first instance based on the 
expressions in Table 1. The nominal value was obtained from the 
known geometric design variables, the solid properties of the 
selected construction materials, and the fluid properties of the 
selected working fluid from the actual NIFTE prototype described in 
detail in Ref. [3[. From this combination of nominal parameters, we 
proceeded to perform a parametric investigation in which each R, C 
and L parameter was varied independently inside a certain range, 
while the rest were kept at their nominal value. The ranges for each 
parameter were selected so as to describe a reasonable incremental 
design perturbation from the existing NIFTE prototype, always 
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leading to an acceptable and realistic configuration that we deem as 
practically achievable in the short term. Table 2 shows the nominal 
values and investigated ranges of each parameter. 

Since all aforementioned components are assumed linear, the 
complete N1FTE network is also linear. For simplicity, but without 
loss of generality, we treat the case of sinusoidal oscillations with 
a single angular frequency to = 2tc f in all quantities. We limit 
ourselves to considering the performance (efficiency) of the NIFTE 
at marginal stability, i.e. in conditions in which it exhibits contin¬ 
uous oscillations of constant amplitude with minimal gain. The 
internal feedback process with gain C, which can be seen on the far 
left in Fig. 2, is increased until this marginal stability condition is 
established. The gain is related directly to the spatial temperature 
gradient established in the source/sink heat exchanger configura¬ 
tion, and thus the difference between the hot and cold tempera¬ 
tures available externally to the device and the design of the heat 
exchanger configuration. In a real NIFTE a gain setting higher than 
that at marginal stability will not result in an unstable system, as 
the system will not be able to sustain continually increasing oscil¬ 
lations and will eventually saturate in a complex non-linear limit 
cycle behavior. Flowever, this is a sub-optimal setting from the 
point of view that it demands a greater temperature difference 
between the heat source and the heat sink (e.g. a higher heat source 
temperature) than that required at the marginal stability condition. 
Thus, the marginal stability set-point can be interpreted as the 
lowest temperature difference condition for which a given NIFTE 
configuration will achieve continuous oscillations. 

At this point we note the relative amplitudes of oscillation and 
(relative) phases of all variables. These are then used in order to 
evaluate the integrals for the heat input rate Q in , the external heat 
supply rate Q h , the net power output W and net exergy flow rateX in 
Eqs. (1—4), which for single-frequency sinusoidal signals simplify to: 


Qi n = (jjT 0 S/n + wTScosX {T,S}/4 

(15) 

Q h = ojT 0 S/n + wf s ScosX {Ts.S }/4 

(16) 

W = P(TcosX{P, U }/2 

(17) 

X = wT s ScosX{Ts,S }/2 

(18) 


where f and f s are the temperature amplitudes of the main cycle 
(i.e. the working fluid) and its thermal environment (measured at 
the internal surface of the source/sink heat exchangers) respec¬ 
tively, S = AS/2 is the entropy amplitude of the main cycle such 

that S = MS is the entropy flow amplitude, and X {a, b} is the phase 
angle between variables a and b. 


Finally, from the quantities in Eqs. (16—18) we can evaluate all 
required efficiencies defined previously. Hence, for the exergetic 
efficiency of the working fluid Vex.wn we obtain the expression, 


_ Wwfl _ Pad COS X {Pad' Qh } 

X - P th cos 4 -{P th ,U th y 


( 19 ) 


while for the exergetic efficiency of the overall device Vex,dev we 
have [21], 


_ Wdev = Pl_ COSXjPT.Ul} 

x p th O t h cos x{Pth,u th } 


( 20 ) 


Similarly, the thermal efficiency of the working fluid ?; th wfl ' s given 

by. 


^th.wfl — 


Wwfl 

Qh 


vjWQjA 
x Vx/ 


^COSXiPad-Uth} 

1 th 


' 2/ir 

,lV7o 


- COS x {Pth - Qfh} 


and the thermal efficiency of the overall device rjth.dev by, 


( 21 ) 


^7th,dev — 


Wdev = 

Qh 

J ^ 

\Ts/T 0 


Wdev/QhV 1 , Pli Ql 

X V x ) p th D th 

l v 1 

+ 2 cosX{P th . Qth) ) . 


cosX{Pi, (Pi} 


( 22 ) 


In the above expressions we have used the results 4.{T s ,S} = 
4{P th , bf t h} an d 4{T,S} = X{P a d. Qh}- Note that even though we 
have made the assumption that the heat transfer processes 
between the heat exchangers and the working fluid, which are 
responsible for driving heat into and out of the cycle, are governed 
exclusively by pure dissipative phenomena in such a way that 
Q°< AT, and even though (consequently) there can be no phase shift 
between S and (T p - T), there is still a phase difference between 
X { T p , S} and X { P, S}. This point and its relation to the performance 
of the NIFTE is investigated closely in Section 3.6. In addition, for 
the two thermal efficiencies r\t h we have used f s /T 0 = 0.05 and 
T s = T 0 , i.e. a 30 K peak difference between the hot source and cold 
sink at T 0 = 300 K. The maximum Carnot efficiency, based on hot 
and cold temperatures of 315 and 285 I< respectively, is 9.5%. 


3. Results and discussion 

This section considers the thermal efficiency of the NIFTE, 
defined as rjth = W/Q h . and the exergetic efficiency of the NIFTE, 


Table 2 

Investigated values of all electrical parameters in Fig. 2 and Table 1. Showing both nominal values that correspond to the actual NIFTE prototype described in detail in Ref. [3]; 
and the full investigated range of each parameter. 


Electrical Parameter 

Symbol and Definition 

Value 

Nominal 

Investigated Range 

Units 

Power cylinder hydrostatic capacitance 

Cp = 

: ndl/Apg 

7.43 x 1(T 8 

3.22 x 10- 9 to 3.22 x 10~ 7 

m 4 s 2 /kg 

Power cylinder liquid mass inductance 

Lp = 

4pl p /7zdp 

3.77 x 10 5 

3.16 x 10 4 to 3.16 x 10 7 

kg/m 4 

Displacer cylinder hydrostatic capacitance 

Q = 

■ Ttd 2 J4pg 

7.35 x 1(T 8 

4.31 x 10~ 8 to 2.58 x 10~ 7 

m 4 s 2 /kg 

Displacer cylinder liquid mass inductance 

id = 

4pl A /-Kd\ 

1.80 x 10 5 

9.96 x 10 4 to 1.39 x 10 6 

kg/m 4 

Feedback connection fluid flow resistance 

R< = 

128 glf/ndf 

2.13 x 10 6 

1.64 x 10 3 to 1.33 x 10 9 

kg/m 4 s 

Feedback connection liquid mass inductance 

L f = 

Aplf/nd^ 

4.74 x 10 6 

4.39 x 10 4 to 3.95 x 10 8 

kg/m 4 

Heat exchanger thermal resistance 

R th = 

= (Sfg/n s ) 2 To/k 

5.02 x 10 s 

2.84 x 10 s to 1.77 x 10 9 

kg/m 4 s 

Adiabatic vapor chamber compressibility 

Qd = 

= ^ao/yPao 

1.76 x 10~ 9 

6.40 x 10-’° to 1.84 x 10~ 8 

m 4 s 2 /kg 

Load fluid flow resistance 

R, = 

128 ph/ndf 

4.08 x 10 6 

4.08 x 10 2 to 5.04 x 10 8 

kg/m 4 s 

Load liquid mass inductance 

L i = 

4pl i /izdf 

1.27 x 10 7 

1.27 x 10 3 to 4.24 x 10 9 

kg/m 4 
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Fig. 3. Effect of power cylinder capacitance C p and inductance L p on the NIFTE device and working fluid thermal efficiencies, T) thdev and Tjth.wfi respectively (left); and device and 
working fluid exergetic efficiencies, i) ex ,dev and i) ex ,wfi respectively (right). The investigated parameters were varied one at a time, with all others set to a nominal value. Abscissas 
normalized by the largest value in the investigated range. 


defined as ?j ex = VV/X, which were evaluated from the expressions 
in Eqs. (19)—(22). As previously stated, efficiencies based on the 
assumption that fluid-mechanical dissipation in the feedback valve 
component is negligible (as in Refs. [1—3]) are referred to as 
working fluid efficiencies ?) w fi and describe the efficiency of the 
thermodynamic cycle undergone by the working fluid; whereas 
efficiencies which account for the fluid-mechanical parasitic 
dissipation in the feedback connection are referred to as device 
efficiencies rjpev and describe the efficiency of the overall device. 
Based on these definitions, results are presented in Figs. 3—8 for the 



effects of selected design parameters (Rs, Cs, Ls in Fig. 2 and Table 1 ) 
on the NIFTE efficiencies: 

(i) Thermal efficiency of the device rjth.devl 

(ii) Thermal efficiency of the working fluid system ?)th,wfi; 

(iii) Exergetic efficiency of the device Tjex.devl and, 

(iv) Exergetic efficiency of the working fluid Vex,wfi- 

In all plots, the parameter under investigation was varied 
independently, with all other parameters set to a nominal value. 



Fig. 4. Effect of displacer cylinder capacitance C d and inductance L d on the NIFTE device and working fluid thermal efficiencies, rj thidev and Tjth.wfi respectively (left); and device and 
working fluid exergetic efficiencies, ?jex,dev and 7j e x.wfi respectively (right). The investigated parameters were varied one at a time, with all others set to a nominal value. Abscissas 
normalized by the largest value in the investigated range. 
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R*, L* [-] R*, L* [-] 

Fig. 5. Effect of feedback valve resistance Rf and connection inductance L; on the N1FTE device and working fluid thermal efficiencies, tjth.dev and ?jth,wfi respectively (left); and device 
and working fluid exergetic efficiencies, ij ex ,dev and ?; eXiW fi respectively (right). The investigated parameters were varied one at a time, with all others set to a nominal value. Abscissas 
normalized by the largest value in the investigated range. 


Specifically, Fig. 3 shows results relating to the role of the power 
cylinder (component 1 in Fig. 2), Fig. 4 shows results relating to the 
displacer cylinder (2 in Fig. 2), Fig. 5 shows results relating to the 
feedback connection (4 and 5 in Fig. 2), and Fig. 6 shows results 
relating to the heat exchangers (10 and 11 in Fig. 2) and the adia¬ 
batic vapor chamber (the total white region 8 in Fig. 2). The 
abscissas (on the horizontal axis) are normalized by the largest 
value in the investigated range. All efficiency values obtained from 
perturbations of the N1FTE from a nominal design can be compared 
to the efficiencies of the nominal NIFTE (i.e. corresponding to the 
actual NIFTE prototype described in Ref. [3]), which are summa¬ 
rized in Table 3. 

By means of comparison, estimated values of rj ex ,dev front 
experiments on a NIFTE water pump [3] with a 125—150 W heater 


sustaining a temperature of 90 °C and a cooling water temperature 
of 4 °C were in the range 0.4—1.7%, while ?)ex,wfi was between 3.7 
and 9.5%. Based on these source and sink temperatures the Carnot 
efficiency is 24%, which corresponds to ?)th,dev = 0.1—0.4% and 
f)th,wfi = 0.9—2.3%. Even though these reported ranges cover 
configurations that extend beyond nominal settings, they provide 
some confidence that the values provided in Table 3 are consistent 
with realistic NIFTE performance characteristics and a valid starting 
point for this current analysis. 

Finally, Figs. 7 and 8 contain results of a sensitivity analysis that 
was performed in order to identify the key parameters that control 
the efficiencies of the NIFTE. Here, the abscissa values of impedance 
on the horizontal axes have been normalized by each variable’s 
nominal value. 




Fig. 6. Effect of thermal resistance R th and adiabatic capacitance C ld on the NIFTE device and working fluid thermal efficiencies, rjth.dev and 7j th wfl respectively (left); and device and 
working fluid exergetic efficiencies, rj ex .dev and t) exw fi respectively (right). The investigated parameters were varied one at a time, with all others set to a nominal value. Abscissas 
normalized by the largest value in the investigated range. 
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Fig. 7. Change in the N1FTE device thermal efficiency >j thds , v due to change in: 
hydrostatic pressure capacitance C p and liquid mass inductance in the power cylinder 
L p ; hydrostatic pressure capacitance Q and liquid mass inductance in the displacer 
cylinder Ly, fluid flow resistance Rf and liquid mass inductance if in the feedback 
connection; thermal resistance in the heat exchangers R th ; and, adiabatic vapor 
chamber capacitance C a d. Abscissa values of impedance Z* on the horizontal axes have 
been normalized by each variable’s nominal value. 

3. J. Power cylinder 

Variations in the thermal efficiency 77 th (left) and exergetic 
efficiency tj e x (right) due to changes in the geometry of the power 
cylinder are shown in Fig. 3. The maximum achievable efficiency 
values obtained by making these changes are tj t h = 5.6% and 
Oex = 76%. These are considerably higher than the values stated in 
Table 3 for nominal NIFTE design and suggest that step improve¬ 
ments are possible. Recall that the capacitance and inductance of 
the power cylinder are given by the expressions C p = nd^/4pg and 
L p = 4p/ p /7td^, respectively. Hence: 

• An increase in the diameter of the power cylinder d p would 
cause C p to increase and L p to decrease; 

• An increase in the liquid density of the working fluid p would 
cause C p to decrease and L p to increase; and, 

• An increase in the time-averaged length of the liquid column in 
the power cylinder l p (related to the time-average of level 6 
in Fig. 2) would cause L p to increase, but would have no effect 
on C p . 

In Fig. 3 we note that the two types of device efficiency (thermal 
and exergetic) lead to consistent design interpretations. For 
example, both p t h And q e x experience a sharp drop as L p is reduced 
below 0.3, and a sharp rise as C p is reduced below 0.2. This simi¬ 
larity can be also observed Figs. 4—6 that relate to the displacer 
cylinder, feedback valve, heat exchangers and adiabatic vapor 
volume. However, the device efficiency is not necessarily optimal 
when the working fluid cycle exhibits its maximum efficiency. For 
example, looking at the figure on the left, rjdev is relatively insen¬ 
sitive to changes in C p except at extremely low values (lower than 
about 0.05) where it shows an improvement, while t; w fl increases 
considerably as C p is gradually reduced below 0.2. This contradic¬ 
tion will be even more evident in Fig. 4 that relates to the effect of 
the displacer cylinder. 

From Fig. 3 it can be seen that C p and L p can be used to make 
significant improvements to the two device efficiencies, pth.dev and 
^ex,dev- In absolute terms, C p can be used to increase tjth.dev by up to 
2.0% and i; e x,dev by up to 26%, while L p can be used to increase tjth.dev 
by up to 4.4% and tjex.dev by up to 61%. From this observation it 



Fig. 8. Change in the NIFTE device exergetic efficiency ij exdev due to change in: 
hydrostatic pressure capacitance C p and liquid mass inductance in the power cylinder 
L p ; hydrostatic pressure capacitance C d and liquid mass inductance in the displacer 
cylinder E d ; fluid flow resistance Rf and liquid mass inductance if in the feedback 
connection; thermal resistance in the heat exchangers R th ; and, adiabatic vapor 
chamber capacitance C, l( |. Abscissa values of impedance Z* on the horizontal axes have 
been normalized by each variable’s nominal value. 

becomes evident that the power cylinder (i.e. the choice of C p and L p ) 
is an important component in determining the overall efficiency of 
the NIFTE. Specifically, the results imply that improved NIFTE effi¬ 
ciencies can be achieved by reducing C p and increasing L p , at least 
within the range of investigate conditions (and given the associated 
assumptions). Both of these improvements can be realized simulta¬ 
neously by decreasing the diameter of the power cylinder d p and 
increasing the liquid density of the working fluid p. An additional 
increase in the length of the liquid column in the power cylinder l p 
would act to further increase L p , without affecting C p . 

Changes to C p and L p can also considerably affect the two 
working fluid efficiencies, pth.wfi and tjex.wfl- Specifically, rjth.wfl 
showed a 5.1% variation over the investigated range of C p , and 
a 5.5% variation over the investigated range of L p . Similarly, over the 
same ranges p e x,wfi varied by 68% and 75%, respectively. 

3.2. Displacer cylinder 

Equivalent plots to those in Fig. 3, but for the displacer cylinder 
are shown in Fig. 4. The capacitance and inductance of the displacer 
cylinder are given by C d = nd^/4pg and L d = 4pl d /v:d^, respec¬ 
tively, such that: 

• An increase in the diameter of the displacer cylinder d d causes 
C d to increase and I d to decrease; 

• An increase in the liquid density of the working fluid p causes 
C d to decrease and L d to increase; and, 

• An increase in the time-averaged length of the liquid column in 
the displacer cylinder / d (related to the time-average of level 7 
in Fig. 2) causes L d to increase, but would have no effect on C d . 

Table 3 

Nominal device and working fluid exergetic and 
thermal efficiencies corresponding to the actual NIFTE 
prototype. 


Efficiency Definition 

Value (%) 

7?th,dev 

0.1 

^th.wfl 

0.3 

V ex.de v 

1.0 

Vex,wfl 

3.9 
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From Fig. 4 it is observed that the magnitude of T^ev is 3—6 times 
smaller than that of r; w fi, or in absolute terms lower by 0.2—4.6%. 
Moreover, the device efficiencies ?jdev are insensitive to the design 
of the displacer cylinder (i.e. the choice of Cd and Ld), with 
^Ith.dev ~ 0.07% and 77ex.dev ~ 0.9% (to within ±1 s.f.) over the entire 
investigated range. On the other hand ?; w n is approximately 
doubled by increasing Cd within the same range. Still, this doubling 
corresponds to a 0.2% absolute increase in ?)th,wfi and an 18% 
absolute increase in rj e x,wfi, which is considerably less than the 
improvement achievable with the careful design of the power 
cylinder suggested previously. 

It is important to recognize that changes to the displacer 
cylinder parameters (Cd and Ld) can have contradictory results on 
the magnitude, but also the trends of the NIFTE’s efficiency, 
depending on the choice of efficiency definition. This highlights 
the fact that the efficiency of the thermodynamic cycle under¬ 
gone by the working fluid ?; w fi can diverge significantly from the 
overall efficiency of the NIFTE device ?7dev. which includes the 
dissipation of useful power in the parasitic feedback valve that is 
external to the main thermodynamic cycle. By extension, it 
demonstrates that neglecting this dissipated power, as was 
assumed in Refs. [1—3] in order to allow a more simple analysis 
and interpretation of the NIFTE, is not generally valid unless the 
feedback resistance Rf is extremely large, in which case the gain 
required in order to attain marginally stable oscillations is also 
very high. In this work we choose ??dev as a more representative 
indication of the overall capabilities of the NIFTE. On this basis, 
we will focus our discussion from this point onwards on the 
overall device efficiency rpi ev , though results will also be pre¬ 
sented for the working fluid efficiency r; w fl. We conclude that the 
displacer cylinder is not important in determining the efficiency 
of the NIFTE. 

3.3. Feedback connection 

Plots concerning the feedback connection are shown in Fig. 5. 
The maximum achievable efficiency values obtained by making 
changes to the feedback connection variables are rjth = 7.1% and 
7]ex = 94%. As with the power cylinder, these values are significantly 
higher than those obtained for the nominal NIFTE design (see 
Table 3), which again implies that the current design can be vastly 
improved with careful design. By comparison, typical thermal effi¬ 
ciencies of around 3—4% were reported for dry free-liquid-piston 
Stirling (Fluidyne) engines [8—11,23,24], though figures as high as 
7% were reported for some larger units. The resistance and induc¬ 
tance of the feedback connection are given by R f = 128 filf/ndf and 
L f = 4plf/itdj, so: 

• An increase in the length of the feedback connection If causes 
both Rf and Lf to increase; 

• An increase in the diameter of the feedback connection df 
causes both Rf and Lf to decrease; 

• An increase in the liquid dynamic viscosity of the working fluid 
p causes Rf to increase, but has not effect on Lf; and, 

• An increase in the liquid density of the working fluid p causes Lf 
to increase, but has no effect on Rf. 

The feedback connection, like the power cylinder, has an 
important role to play in determining the eventual efficiency of the 
NIFTE. An examination of the effects of the parameters Rf and Lf on 
our selected efficiency definition Tjdev (i.e. only the circles) in Fig. 5, 
reveals that changes to Rf have a much greater ability to influence 
Pdev than changes to Lf. In particular, Fig. 5 (left) shows that Rf (filled 
circles) can be used to improve jjth.dev by up to 3.9% in absolute 
terms, and (right) to improve 7/ e x,dev by up to 51% in absolute terms. 


The corresponding values for Lf (hollow circles) are 0.6% and 9%, 
respectively. Thus, the sensitivities of 7) t h,dev and pex.dev to Rf are 
about 6 times greater than their sensitivities to Lf. The sensitivities 
of the important device efficiency definitions t] d ev to all NIFTE 
parameters will be examined in detail in the discussion relating to 
Figs. 7 and 8. 

So the model predicts that an efficiency improvement can be 
achieved by reducing Rf, and to a lesser extent Lf. Both of these can 
be realized simultaneously by decreasing the length If and 
increasing the diameter df of the feedback connection. An addi¬ 
tional decrease in the liquid dynamic viscosity of the working fluid 
p and a decrease in the liquid density of the working fluid p will also 
act to independently lessen Rf and Lf, respectively. 

Recall however, from the sub-section on the power cylinder 
optimization, that an increased liquid density p was required there 
for maximum efficiency. Hence, the optimization of the NIFTE with 
respect to this variable is more complex. One approach to avoid this 
conflict in requirements is to set Rf and Lf by designing the geom¬ 
etry of the feedback connection (If and df), and to set C p and L p 
independently by designing the geometry of the power cylinder 
(d p ), without the need to change the working fluid. 

3.4. Heat exchangers and adiabatic vapor chamber 

Finally, we turn to the heat exchangers and the total vapor 
volume at the top of the NIFTE. Plots concerning the effect of these 
two components on the efficiency of the NIFTE are shown in Fig. 6. 
The heat transfer process in the heat exchangers is described by 
the thermal resistance R th = (s fg /vf g ) 2 T 0 /k, while the adiabatic 
capacitance associated with the adiabatic vapor chamber is 
C a d = Vao/yPzo- From these expressions, an increase in R t h can be 
established with either a: 

• Higher specific enthalpy due to phase change Sf g ; or, 

• Lower specific volume due to phase change Vf g ; or, 

• Higher time-averaged (saturation) temperature T 0 ; or, 

• Lower thermal conductance k. 

In addition, C a d can be increased with either a: 

• Higher mean total volume in the vapor phase V ao ; or, 

• Higher ratio of specific heat capacities of the working fluid in 
the vapor phase y; or, 

• Lower time-averaged (saturation) vapor pressure P ao , which is 
directly linked to the (saturation) temperature T 0 . 

As before, we examine the effects of the two independent 
variables (here R t h and C a d) on our selected (device) efficiency 
definition Tjdev (he. only the circles) in Fig. 6. The thermal resistance 
Rth (filled circles) has a less noticeable effect on ?jdev than that of C a d 
(hollow circles). In fact, in Fig. 6 the sensitivities of both ??th,dev and 
Vex.dev to C a d are approximately 5 times greater than their sensi¬ 
tivities to R t h, so the first priority in maximizing rj d ev should be to 
maximize C a d. In more detail, decreasing R t h within the investigated 
range increases t| t h,dev by 0.1% (or a factor of about 2), and 
increasing C a d increases ?7th,dev by 0.5% (or a factor close to 5). 
Similarly, decreasing R t h increases rj e x.dev by 1.4%, and increasing C a d 
increases 7) e x,dev by 6.9%. 

The selection of the working fluid is an important, but compli¬ 
cated issue in itself and depends on a number of additional factors 
(not contained in the current study) that must be weighed before 
a decision can be made. Such factors include for example its boiling 
point in relation the hot Th and cold T c temperatures, amongst 
others. Specifically with respect to the boiling point, one basic 
requirement is that this must be somewhere between Th and T c if 
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the working fluid is going to evaporate and condense, otherwise the 
NIFTE will not be able to operate. 

Given that there is little scope in affecting the specific heat 
capacity ratio of the working fluid in the vapor phase y, it is only 
possible to establish considerably larger values of C a d by increasing 
the total volume in the vapor phase V ao , and/or decreasing the time- 
averaged (saturation) pressure P ao . With respect to the former, the 
total adiabatic vapor chamber volume V ao can be increased by using 
power and displacer cylinders with larger diameters and lengths, as 
well as vapor connection pipes with larger diameters and lengths. 
Yet, the requirement concerning the power cylinder length is in 
partial contradiction with the condition that was established in the 
power cylinder sub-section for maximum efficiency, i.e. that 
a narrow (small d p ) power cylinder ought to be preferred. One way 
to overcome this problem is to use a long, but small diameter power 
cylinder, thus achieving both small C p and large C a( j. 

In contrast to V ao , the designer has little control over the time- 
averaged (saturation) pressure P ao . This pressure appears as an 
external boundary condition observed by the NIFTE, as it is deter¬ 
mined by the time-mean pressure in the load to which the NIFTE is 
connected. In the linear model employed here the time-averaged 
temperature T 0 that is found half-way between the source and 
sink temperatures, which also coincides with saturation conditions. 
Since P ao is coupled to T 0 , the time-averaged pressure P ao also 
determines the temperatures in the heat exchangers. In fact, we 
expect that lower P ao will be associated with lower heat exchanger 
temperatures. Nevertheless, we expect that C ac | will be higher 
under these conditions, leading to improved efficiencies. A lower T 0 
will also help to decrease R t h that can have a secondary beneficial 
effect on the NIFTE efficiency. Note however, that the linear analysis 
performed here is too simplistic to capture non-linear effects that 
become increasingly important at higher P ao and higher heat 
exchanger temperatures. 

3.5. Sensitivity analysis 

In this section we re-examine the data presented previously, but 
in the form of a sensitivity analysis in order to identify the R, C and L 
parameters (and NIFTE variables) that most strongly affect the 
efficiencies of the NIFTE. We focus here only on the overall device 
efficiencies jjdev. which include parasitic fluid power dissipation 
due to viscous and pressure drag in the feedback connection (due to 
R t h). Throughout this section we assume that the designer has no 
control over the load which the device interfaces with, as this is set 
externally by the application. We therefore set the load to its 
nominal value and exclude all load parameters from consideration. 


Fig. 7 (with key values summarized in Table 4) contains plots of 
the changes in the thermal efficiency of the device rjth.dev due to 
changes in all independent device parameters, including: the 
capacitance due to hydrostatic pressure C p and inductance due to 
liquid mass (inertia) in the power cylinder L p ; the capacitance due 
to hydrostatic pressure Q| and inductance due to liquid mass 
(inertia) in the displacer cylinder Ld; the resistance due to fluid flow 
(drag) Rf and inductance due to liquid mass (inertia) Lf in the 
feedback connection; the thermal resistance in the heat exchangers 
Rthi and, the capacitance associated with compressibility of the 
adiabatic vapor chamber C a d. 

Fig. 8 (with related key results summarized in Table 5) is iden¬ 
tical to Fig. 7, but focuses on changes to the exergetic efficiency of 
the device i) e x.dev- It is possible to conclude from Figs. 7 and 8, as 
well as the accompanying Tables 4 and 5, that the efficiencies of 
NIFTE are most sensitive to changes in the feedback connection 
(both Rf and Lf), followed by the power cylinder (C p and I p ), and to 
a lesser extent the thermal resistance in the heat exchangers R t h 
and the capacitance of the adiabatic vapor chamber C ai j. Based on 
this information it becomes evident that the selection of the 
feedback valve as the natural operating point tuning component is 
the right approach. 

3.6. Temperature-entropy diagrams 

It has been shown in the previous sections that the feedback 
connection and power cylinder are important components, whose 
optimization can lead to significant improvements in efficiency 
with respect to an existing (nominal) NIFTE prototype. For instance, 
a decrease in Rf resulted in an improvement in rjth.dev from 
a nominal 0.1% up to 3.9% in absolute terms, and similarly in ??ex,dev 
from a nominal 1% up to 52%. In another example, a decrease in C p 
was capable of increasing ?)th,dev from a nominal 0.1% to 2.1%, and 
Vex,dev from a nominal 1% to 27%. Corresponding working fluid 
efficiencies also showed considerable improvements, with Tjth.wfl 
increasing from 0.3% to 7.1% at the optimal Rf and 5.2% at the 
optimal C p , respectively; and with rjex.wfl increasing from 4% to 94% 
at the optimal Rf and 69% at the optimal C p setting. 

In order to gain some insight into the underlying reasons for 
these vast potential improvements, it is useful to examine the main 
(internally reversible) thermodynamic cycle undergone by the 
working fluid on a T—S diagram and also to compare this with the 
greater area that represents the full (thermally irreversible) cycle as 
experienced at the external heat source and sink. As explained 
previously, the distinction between these two rests on the choice 
of temperature, either T or T s , and the area between the two 


Table 4 

Device thermal efficiency 7j t h,dev sensitivity coefficients corresponding to Fig. 7. Showing both nominal coefficients that correspond to the actual NIFTE prototype described in 
detail in Ref. [3]; and the maximum (of the modulus) of each coefficient over the full investigated range of each parameter, as presented in Table 2. 


Parameter 

Symbol and Definition 

Device Efficiency Sensitivity Coefficient 

Nominal Maximum (Modulus) 

Power cylinder hydrostatic capacitance 

Cp = 

= Kdl/4pg 

2.92 

X 

10" 3 

1.85 

X 

io 1 n 

Power cylinder liquid mass inductance 

Ip = 

4pl P /Kdl 

3.05 

X 

1(T 4 

2.26 

X 

lO" 3 

Displacer cylinder hydrostatic capacitance 

Q = 

= ndl/4pg 

4.06 

X 

1<T 4 

4.21 

X 

10- 4 

Displacer cylinder liquid mass inductance 

Ld = 

= 4pk/t tdjj 

2.66 

X 

10" 5 

4.00 

X 

10" 5 

Feedback connection fluid flow resistance 

Rf = 

= 128 filf/ird/ 

1.70 

X 

10" 3 

1.45 

X 

io +1 [- 1 ] 

Feedback connection liquid mass inductance 

L { = 

; 4pl f /Kdj 

7.41 

X 

1(T 4 

1.45 

X 

io- 1 [ b ] 

Heat exchanger thermal resistance 

^th 

= (Sfg/yfg) T 0 /k 

2.63 

X 

10“ 3 

3.46 

X 

10" 3 

Adiabatic vapor chamber compressibility 

Cad 

= Vao/yPao 

9.65 

X 

10- 4 

3.32 

X 

IO" 3 


a Primary control component major variable. 
b Primary control component minor variable. 
c Secondary control component main variable. 
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Table 5 

Device exergetic efficiency ij ex .dev sensitivity coefficients corresponding to Fig. 8. Showing both nominal coefficients that correspond to the actual NIFTE prototype described in 
detail in Ref. [3]; and the maximum (of the modulus) of each coefficient over the full investigated range of each parameter, as presented in Table 2. 


Parameter 

Symbol and Definition 

Sensitivity Coefficient 

Nominal 

Maximum 

Power cylinder hydrostatic capacitance 

Cp = 

S2 

(NO. 

£ 

1.19 

X 

10" 3 

1.91 

X 

io° n 

Power cylinder liquid mass inductance 

Lp = 

= 4p/ p /7td 2 

1.05 

X 

10- 4 

2.79 

X 

10" 2 

Displacer cylinder hydrostatic capacitance 

Cd = 

= rafji/4 Pg 

1.46 

X 

10“ 3 

3.41 

X 

10" 3 

Displacer cylinder liquid mass inductance 

id = 

= 4pl d /izd 2 

1.68 

X 

10- 4 

2.47 

X 

10- 4 

Feedback connection fluid flow resistance 

Rf = 

- 128 fdf/TtdJ! 

2.02 

X 

1(T 4 

1.93 

X 

io+ 2 [ a ] 

Feedback connection liquid mass inductance 

L f = 

= 4p!f/mfJ 

3.97 

X 

10" 3 

1.66 

X 

10° [ b ] 

Heat exchanger thermal resistance 

^th 

= (Sfg/vfg) T 0 /k 

8.89 

X 

10" 3 

2.38 

X 

10" 2 

Adiabatic vapor chamber compressibility 

Q d 

— Vao/yPao 

1.97 

X 

io- 3 

9.26 

X 

10- 3 


a Primary control component major variable. 
b Primary control component minor variable. 
c Secondary control component main variable. 


represents a power loss (exergy destruction) through irreversible 
heat transfer that occurs in R t h between the source/sink and the 
working fluid. Fig. 9 shows results from this investigation for the 
effects of Rf (left) and C p (right) on the T—S (solid lines, filled 
squares) and T s —S (dotted lines, hollow squares) NIFTE diagrams 
undergone in the relative to the nominal configuration (circles). 

The areas within the two ovals are a direct indication of the 
efficiencies of the corresponding cycle. As expected, it can be 
observed in Fig. 9 that the T s —S diagrams (hollow symbols) always 
enclose the T—S diagrams (filled symbols). It can also be seen that 
the T s —S diagrams cover the whole vertical range (height) in the 
graphs T' that corresponds to the amplitude of temperature oscil¬ 
lations in the heat exchangers f s = 15 K, again as expected. The 
nominal design however has a greatly reduced working fluid 
temperature oscillation amplitude t = 1 K, revealing a substantial 
exergy loss in the thermal resistance R t h. This helps to explain the 
relatively low ?jex,wfi and ?)ex,dev. since both of these efficiency 
measures are adversely affected by these exergetic losses. 

Continuing, the most important observation that can be made 
based on the plots in Fig. 9 concerns the effects on the T-S diagrams 
(solid lines, square symbols) of optimizing either the value of Rf 
(left) or C p (right). Compared to the nominal NIFTE, the optimized 
designs are associated with much larger amplitudes in the working 
fluid temperature oscillations. Specifically, we can see clearly from 
these plots that the underlying effect of the parametric optimiza¬ 
tion process is to increase T to 141< when Rf is optimized and to 10 K 
when C p is optimized. This corresponds directly to the higher rj ex in 
the optimal designs relative to the nominal NIFTE. In addition, the 
larger relative increase in f that can be generated by varying Rf 



rather than C p within their perturbation envelopes can explain the 
higher i) ex of the Rf optimal NIFTE design. These values are 
summarized in Table 6. 

As a final comment, it can be observed in Table 6 that unlike the 
nominal NIFTE, the optimized designs are capable of keeping the 
phase between 4.{T s ,S} and 4{T, S} low (see Eqs. (19)—(22)). This 
phasing is an important parameter in determining the eventual 
efficiency as it controls (together with the amplitudes f s and T) the 
area of the Lissajous ovals in Fig. 9, with zero phase corresponding 
to maximum efficiency and 90° phase corresponding to zero effi¬ 
ciency. However, close inspection of various other sub-optimal 
designs (not shown here) reveals even lower phases than those 
reported in Table 6, suggesting that the primary reason for effi¬ 
ciency improvement at lower Rf and C p is the increase in temper¬ 
ature amplitude ratios, rather than a decrease in the phase 
difference between 4{T s ,S} and 4{T,S}. 


Table 6 

Result of optimization of the NIFTE with respect to (w.r.t.) Rf and C p in terms of 
external and working fluid temperature amplitudes (T s and T) and corresponding 
phases ( 4 {T s ,S} and 4 {T,S}). Comparison with the nominal design. 


Variable 

Nominal Value 

Optimized NIFTE 

w.r.t. Rf 

w.r.t. C p 

Ts (K) 

15 

15 

15 

T (K) 

1 

14 

10 

4 {T s ,S} (deg.) 

3° 

8° 

4° 

4 {T,S} (deg.) 

52° 

8° 

6° 



Fig. 9. First-order linear approximation phase plane plots (Lissajous ovals) of the NIFTE thermodynamic cycle, showing the nominal configuration (circles) with 
Rf = 2.13 x 10 6 kg m 4 s 1 and C p = 7.43 x 10 3 m 4 s 2 kg~\ as well as optimized cycles with respect to R f (=1.64 x 10 3 kg m 4 s ') with C p fixed at its nominal value (left, square 
symbols) and C p (=3.22 x 10 9 m 4 s 2 kg -1 ) with Rf fixed at its nominal value (right, square symbols). Each figure includes both working fluid T-S (solid lines, filled symbols) and 
device T s -S (dotted lines, hollow symbols) diagrams. Abscissa values of entropy fluctuations S t f,* on the horizontal axes are normalized by their amplitude. 
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4. Conclusions 

Starting from a first-order linear dynamic model of the N1FTE 
that consists of a network of interconnected spatially lumped 
components, we investigated parametrically the effect of various 
component variables (geometric and other) on the thermal (first 
law) and exergetic (second law) efficiencies of this device. Results 
from the parametric optimization process in which the parameters 
were varied independently within a specified range, were then 
compared to efficiency values associated with a nominal N1FTE 
configuration, based on values that were estimated from an exist¬ 
ing prototype. Great improvements were demonstrated relative to 
the nominal NIFTE design, implying that the careful design can 
have a significant effect on the performance of the final device. 

The study highlighted critical components that require careful 
optimization, namely the feedback connection valve, the power 
cylinder, the adiabatic volume and the thermal resistance in the 
heat exchangers. Relative to the nominal NIFTE, an efficient design 
would feature a lower feedback connection valve resistance, with 
a shorter connection length and larger connection diameter; 
a smaller diameter but taller power cylinder; a larger (time-mean) 
combined vapor volume at the top part of the device; and improved 
heat transfer behavior (i.e. reduced thermal resistance) in the hot 
and cold heat exchanger blocks. 

The effect of these modifications was to increase the exergetic 
efficiency of the device by 50% points, which also corresponds to 
a 3.8% point increase in thermal efficiency. The NIFTE was modeled 
as operating across a temperature difference between heat source 
and sink that one would associate with low grade heat (301<). Given 
a maximum theoretical (Carnot) efficiency associated with these 
temperatures of 9.5%, the reported improvement is considered 
noteworthy. Close inspection of T—S diagrams revealed that the 
underlying reason for the improvement was an increase in the 
working fluid temperature amplitude relative to that made avail¬ 
able to it by the heat source and sink, as it undergoes the ther¬ 
modynamic cycle. 

Acknowledgments 

The authors would like to thank Mr. Y. Malik and Miss R. Solanki 
for their assistance in generating the data. 

References 

[1] Smith TCB. Power dense thermofluidic oscillators for high load applications. 
AIAA-2004-5758. In: Energy conversion engineering. IECEC 2004: Proceedings 
of the 2nd International energy conversion Engineering Conference; 2004 Aug 
16-19; Providence (RI). USA. Reston (VA): AIAA; 2004. p. 1-15. 


[2] Smith TCB. Asymmetric heat transfer in vapour cycle liquid-piston engines. In: 
Stirling engine. ISEC 2005: Proceedings of the 12th International Stirling 
engine Conference and technology Exhibition; 2005 Sept 7—9. Durham, UK. 
Durham: Durham University; 2005. p. 302-14. 

[3] Smith TCB. Thermally driven oscillations in dynamic applications, PhD thesis. 
Cambridge: University of Cambridge, UK; 2006. 

[4] Backhaus S, Swift GW. A thermoacoustic-stirling heat engine: detailed study. 
J Acoust Soc Am 2000;107:3148-66. 

[5] Wheatley J, Hofler T, Swift GW, Migliori A. An intrinsically irreversible 
thermoacoustic heat engine. J Acoust Soc Am 1983;74:153-70. 

[6] Wheatley J, Hofler T, Swift GW, Migliori A. Experiments with an intrinsically 
irreversible acoustic heat engine. Phys Rev Lett 1983;50:499-502. 

[7] Ceperley PH. A pistonless Stirling engine—The travelling wave heat engine. 
J Acoust Soc Am 1979;66:1508-13. 

[8] West C. The Fluidyne heat engine. Research report. Harwell: atomic energy 
Research Establishment; 1971. Report No. AERE-R 6775. 

[9] West CD, Pandey RB. Laboratory prototype Fluidyne water pump. In: Energy 
conversion engineering. IECEC 1981: Proceedings of the 16th Intersociety 
energy conversion Engineering Conference; 1981 Aug 9—14; Atlanta (GA). 
USA. New York (NY): ASME; 1981. p. 1916-8. 

[10] West CD. Dynamic analysis of the Fluidyne. In: Energy conversion engi¬ 
neering. IECEC 1983: Proceedings of the 18th Intersociety energy conversion 
Engineering Conference; 1983 Aug 21-26; Orlando (FL). USA. New York (NY): 
AIChE; 1983. p. 779-84. 

[11] Stammers CW. The operation of the Fluidyne heat engine at low differential 
temperatures. J Sound Vib 1979;63:507-16. 

[12] Redlich RW, Berchowitz DM. Linear dynamics of free-piston Stirling engines. 
Proc Inst Mech Eng Part A Power Process Eng 1985;199:203-13. 

[13] Walker G, Senft JR. Free piston Stirling engines. Berlin: Springer; 1985. 

[14] Wood JG, Lane NW. Advanced 35 We Stirling engine for space power appli¬ 
cations. In: El-Genk MS, editor. STAIF 2003: Proceedings of the space tech¬ 
nology and applications International Forum; 2003 Feb 2—5; Albuquerque 
(NM). USA. Melville (NY): AIP; 2003. p. 662-7. 

[15] Huang BJ, Chuang MD. System design of orifice pulse-tube refrigerator using 
linear flow network analysis. Cryog 1996;36:889-902. 

[16] Payne PR, Brown RG, Brown JP. Water pulsejet research. Final report. 
Annapolis (MD): Payne Inc; 1979 Apr. Report No. DOE/ET/20288—Tl. Contract 
No.: AC01-77ET20288. Sponsored by the Department of Energy. 

[17] Kentfield JAC. Fundamentals of idealized airbreathing pulse-detonation 
engines. J Propuls Power 2002;18:78-83. 

[18] Kentfield JAC. Thermodynamics of airbreathing pulse-detonation engines. 
J Propuls Power 2002;18:1170-5. 

[19] Organ AJ. Stirling and pulse-tube cryo-coolers. 1st ed. Bury St. Edmunds: John 
Wiley and Sons; 2005. 

[20] Shukuya M, Hammache A. Introduction to the concept of exergy — for a better 
understanding of low-temperature-heating and high-temperature-cooling 
systems. Report. Espoo: VTT Technical Research Centre of Finland; 2002 
Apr. Research Notes 2158. 

[21 ] Solanki R, Galindo A, Markides CN. Dynamic modelling of a two-phase ther¬ 
mofluidic oscillator for efficient low grade heat utilization: effect of fluid 
inertia. J Appl Energy 2012;89:156-63. 

[22] Unsal B, Ray S, Durst F, Ertung O. Pulsating laminar pipe flows with sinusoidal 
mass flux variations. Fluid Dyn Res 2005;37:317-33. 

[23] Goldberg LF, Rallis CJ. A prototype liquid-piston free-displacer Stirling engine. 
In: Energy conversion engineering. IECEC 1979: Proceedings of the 14th 
Intersociety energy conversion Engineering Conference; 1979 Aug 5—10; 
Boston (MA). USA. Washington (DC): ACS; 1979. p. 1103-8. 

[24] Elston MJ, Lurie MS, Rallis CJ, Kilgour DB. Further development of the 
Fluidyne liquid-piston engine. In: Energy conversion engineering. IECEC 
1982: Proceedings of the 17th Intersociety energy conversion Engineering 
Conference; 1982 Aug 8—12; Los Angeles (CA). USA. New York (NY): IEEE; 
1982. p. 1755-60. 




